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Abstract.  The  Method  of  Dimensionality  Reduction  (MDR)  is  a  method  of  calculation 
and  simulation  of  contacts  of  elastic  and  viscoelastic  bodies.  It  consists  essentially  of 
two  simple  steps:  (a)  substitution  of  the  three-dimensional  continuum  by  a  uniquely 
defined  one-dimensional  linearly  elastic  or  viscoelastic  foundation  (Winkler  foundation) 
and  (b)  transformation  of  the  three-dimensional  profile  of  the  contacting  bodies  by 
means  of  the  MDR-transformation.  As  soon  as  these  two  steps  are  completed,  the 
contact  problem  can  be  considered  to  be  solved.  For  axial  symmetric  contacts,  only  a 
small  calculation  by  hand  is  required  which  does  not  exceed  elementary  calculus  and 
will  not  be  a  barrier  for  any  practically -oriented  engineer.  Alternatively,  the  MDR  can 
be  implemented  numerically,  which  is  almost  trivial  due  to  the  independence  of  the 
foundation  elements.  In  spite  of  their  simplicity,  all  the  results  are  exact.  The  present 
paper  is  a  short  practical  guide  to  the  MDR. 
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1.  Introduction 

In  the  recently  published  book  [1],  the  so-called  method  of  dimensionality  reduction 
(MDR)  is  described  for  the  first  time  in  detail.  MDR  can  be  traced  back  to  the  solution  of 
the  normal  contact  problem  by  Galin  (Russian  Academy  of  Sciences)  in  the  1940s  [2]. 
His  results  were  later  published  by  Sneddon  and,  in  this  way,  made  public  to  the  western 
world  [3].  The  method  of  dimensionality  reduction  takes  these  results  and  puts  them  into 
such  a  form  that  even  a  layman  in  the  field  of  contact  mechanics  can  use  them  for  a 
multitude  of  contact  mechanical  problems.  In  doing  this,  it  merges  the  ideas  and  results 
from  Cattaneo  [4],  Mindlin  [5],  Jaeger  [6],  and  Ciavarella  [7]  about  a  close  relationship 
between  normal  and  tangential  contacts,  the  solutions  of  Galanov  and  Borodich  [8,  9,  10] 
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for  adhesive  contacts  of  axially-symmetric  profiles  of  power  functions  (later  found 
independently  by  Yao  and  Gao  [11]),  as  well  as  the  theory  of  Lee  and  Radok  about  the 
relationship  between  elastic  and  viscoelastic  contacts  [12,  13]. 

The  book  [1]  contains  all  of  the  necessary  evidence  and  many  examples  of  how  to  apply 
the  MDR.  However,  it  has  proven  to  be  too  comprehensive  for  practical  users.  There  is  a 
need  for  the  fundamental  ideas  and  "recipes"  of  the  MDR  to  be  presented  in  a  concise  way 
without  extensive  reasoning  or  proof,  a  sort  of  "user's  handbook."  This  work  is  dedicated  to 
exactly  such  a  practical  instruction  for  the  method  of  dimensionality  reduction. 


2.  Two  Introductory  Steps  of  the  MDR 

We  consider  a  contact  between  two  elastic  bodies  with  moduli  of  elasticity  of  Ej  and  E2, 
Poison's  numbers  of  Vj  and  v2,  and  shear  moduli  of  Gj  and  G2,  accordingly.  In  this  work,  we 
restrict  ourselves  to  the  axially-symmetric  profiles,  which  is  not  necessarily  required.  A 
generalization  to  profiles  that  are  not  axially-symmetric  is  possible,  but  is  not  considered  in 
this  work.  We  denote  the  difference  between  the  profiles  of  bodies  as  z  =  f(r).  In  the 
framework  of  the  MDR,  two  independent  steps  are  conducted: 

The  first  step:  First,  the  three-dimensional  elastic  (or  viscoelastic)  bodies  are  replaced 
by  a  one-dimensional  linearly  elastic  foundation.  This  is  considered  to  be  a  linear  array  of 
elements  having  independent  degrees  of  freedom  and  a  sufficiently  small  separation 
distance  Ax,  Fig.  1 . 


XXX 


a)  b)  c) 

Fig.  1  One-dimensional  foundation  of  different  materials:  elastic  foundation  (a),  purely 
viscous  foundation  (b),  and  viscoelastic  foundation  (c)  with  an  example  rheology 
according  to  Kelvin- Voigt 


In  the  simplest  case  of  the  elastic  contact,  the  foundation  consists  of  linearly  elastic 
spring  elements  that  have  normal  stiffness  Akz  and  tangential  stiffness  Akx  (Fig.  la): 

1      1  _    2    1  —  2 
Akz  =  E*Ax  with  —  =  -^i-  +  -^-2- ,  (1) 

M,  =  G" Ax  with  -4  =  7—^~  +  7—^L  ■  (2) 
G       AG,  4G2 

Starred  values  E*  and  G*  denote  the  effective  elastic  moduli.  Incompressible  linearly 
viscous  materials  are  presented  by  a  linear  damping  element  with  damping  coefficient  Ay 
(Fig.  lb),  which  is  dependent  on  the  viscosity  t]  of  the  viscoelastic  partner  according  to: 
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A/  =  4rjAx . 


(3) 


Arbitrary  combinations  of  these  two  base  elements  are  also  possible  in  order  to  satisfy 
the  most  complicated  elastomers  -  Fig.  lc,  for  example,  shows  a  viscoelastic  foundation 
built  out  of  elements  of  in  parallel  connected  springs  and  dampers  (Kelvin- Voigt  model). 
In  this  paper,  we  will  restrict  ourselves  to  the  case  of  "elastically  similar"  materials: 


1  — 2v,     1  —  2v2 


(4) 


which  guarantees  the  independence  of  the  normal  and  tangential  contact  problems  [14]. 
This  condition  is  always  met  in  important  cases  of  contacts  between  the  bodies  with  the 
same  elastic  properties  or  those  between  a  rigid  body  and  an  elastomer. 

The  second  step:  In  the  second  step,  three-dimensional  profile  z  =fC/)  (Fig.  2,  left)  is 
transformed  into  a  one-dimensional  profile  (Fig.  2,  right)  according  to: 


g(x)  =  |*|  J 


fir) 


o 


rdr. 


The  reverse  transformation  is: 


f(r)  =  -l 

IT  •> 


2'f  g(x) 


no4 


2  2 

r  —x 


-dx . 


(5) 


(6) 


f(r)  ^> 


Fig.  2  The  three-dimensional  profile  is  transformed  into  a  one-dimensional  profile  using  the  MDR 


For  a  less  trivial  example,  we  consider  the  contact  of  a  parabolic  profile  with  a  worn 
tip  (Fig.  3): 
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r  —b 
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for  0<r<b 
for  b<  r  <a 


(7) 


The  MDR  transformed  profile  according  to  Eq.  (5)  is  given  by: 

0  for  0<\x\<b 

gix)=* 


R 


- -Jx2  —  b2     for  b  <  \x\  <  a 


(8) 
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§U,9) 


original  profile  (3D) 
equiviLleiit  profile  ( ID) 


Fig.  3  Parabolic  indenter  with  a  "worn"  tip:  original  (solid  line)  and  equivalent 
(dashed  line)  profiles  for  comparison 


Examples  for  the  MDR  transformation 

By  inserting  profiles  into  Eq.  (5)  that  correspond  to  a  cylinder,  paraboloid,  cone,  or  an 
arbitrary  power  function  z  ~  r'\  we  obtain  the  MDR  transformed  one-dimensional  profiles 
which  are  summarized  in  Table  1,  where: 


k.  = 


r  nT\n 
VJi       1 2 


— - — ^-  and  T(n)  :=  \t"  [e  'dt  is  the  gamma  function, 
nil  n 


r  -  + 

2  2 


Table  1  The  three-dimensional  profile  is  transformed  into  a  one -dimensional  profile  using 
the  MDR 


^  j  arbitrary 

.J  n 


0,  r<a 
co,  r  >  a 
0,  x  <  a 
oo,  x  >  a 


r2  /2R 
x2  (R 


rtanf? 


—  dtanf?  Kn^?J■'^, 
2  I  I  1 


Eq.  (7) 


Eq.  (8) 


3.  Calculation  Step  of  the  MDR  Using  the  Example  of 
a  Normal  Contact  without  Adhesion 

The  one-dimensional  profile  according  to  Eq.  (5)  is  now  pressed  into  an  elastic 
foundation  corresponding  to  Eq.  (1)  with  normal  force  FN  (see  Fig.  4).  The  normal 
surface  displacement  at  point  x  within  the  contact  area  results  from  the  difference  between 
indentation  depth  d  and  profile  form  g: 
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Fig.  4  Equivalent  model  for  a  Hertzian  contact 

uz(x)  =  d-g(x).  (9) 

At  the  edge  of  non-adhesive  contact  x  =  +  a,  the  surface  displacement  must  be  zero: 

u_(±a)  =  0  =>  d:=g(a).  (10) 

This  equation  determines  the  relationship  between  the  indentation  depth  and  contact 
radius  a.  It  should  be  noted  that  this  relationship  is  independent  of  the  rheology  of  the 
medium.  The  force  of  a  spring  at  the  point  x  is  proportional  to  the  displacement  at  this  point: 


AF  (x)  =  Ak7u7  (x)  -  E  u  (x)Ax , 


(11) 


and  the  sum  of  all  spring  forces  must  correspond  to  the  normal  force  in  equilibrium.  In  the 
limiting  case  of  very  small  spring  separation  distances  Ax  — >  dx,  the  summation  becomes 
the  integral: 


FN  :=  E"  J  uz(x)dx  =  2E  j"(c? -g(x))dx  . 


(12) 


Equation  (12)  provides  the  normal  force  in  dependence  on  the  contact  radius  and  on  the 
indentation  depth,  if  Eq.  (10)  is  taken  into  account. 
We  now  define  linear  force  density  qz(x): 


AFJx) 

q ,  (x)  :=  =  E  u,  (x)  . 

Ax 


(13) 


The  stress  distribution  in  the  original  three-dimensional  system  can  be  determined 
with  the  help  of  the  one-dimensional  distributed  load  using  the  integral  transformation 
[15]: 

q(x) 


1 

P(r)  =  —  | 

77"  J 


2  2 

x  —  r 


-  dx  . 


(14) 


The  normal  surface  displacement  (both  inside  and  outside  of  the  contact  area)  is  given  by 
the  transformation: 


u  Ax) 


-  Ax . 


(15) 
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We  point  out  once  more  that  one-dimensional  values  contain  the  independent  variable 
x,  while  three-dimensional  values  contain  radial  variable  r;  function  uz(x)  inside  the 
integrand  of  Eq.  (15)  means  the  surface  displacement  of  the  linearly  elastic  foundation. 
We  will  carry  out  the  last  two  transformations  using  the  example  of  a  cone.  The  three- 
dimensional  profile  in  this  case  is  f(r)  =  ntan&,  the  MDR  transformed  profile  is 


n 


n  i 


g(x)  =  —\x\tan0 .  The  displacement  within  the  contact  area  is  u  (x)  =  d  be  tan#.  The 


linear  force  density  is  q(x)  =  E"  I  d  —  —  \x\taa0  I  and  its  derivative  is  q'(x)  =  --E*  tan# 


n 


(16) 


(for  positive  x).  Insertion  into  Eq.  (14)  and  Eq.  (15)  results  in: 


p(r)  =  —  tan«9  \-j== 

Z  r  VX  - 

2  atd-{nl2)x\zn0 
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=  =  —  tannin 
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2d 
ax  =  — 
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u 
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b: 

a  \ 

j 

(17) 


Equation  (17)  gives  the  normal  surface  displacement  outside  of  the  contact  area.  A 
similar  calculation  for  a  parabolic  profile  f(r)  =  r2/(2R)  initially  provides  uz(x)  =  d-x2/R 
and,  after  insertion  into  Eq.  (15),  results  in: 

d 


uAr)-. 


71 


'  r 

i; 

• arcsin  j 

V 

(18) 


For  the  contact  of  a  flattened  paraboloid,  Eq.  (7),  we  obtain  the  contact  radius  by 
using  Eq.  (10): 


«,(a)  =  0  =>  d  :=  g(a)  =  —4a2  -b2  . 

R 


and  the  normal  force  using  Eq.(12): 

2E"  r 
R 

0  JV  0 


FN  =2E'jddx  jx^lx2-b2dx  =  (2a2  +b2)  ■  s/a2 -b 


(19) 


(20) 


These  results  correspond  of  course  with  those  of  Eijke  [18]  from  the  three-dimensional 
theory. 

Examples  for  normal  contacts 

Insertion  of  the  MDR  transformed  profiles  into  Eqs.  (10)  and  (12)  and  an  elementary 
integration  provides  the  results  for  the  "classical  profiles"  of  a  cylinder  [17],  sphere  [16], 
and  cone  [3]  as  well  as  the  general  power  profile  [2],  which  are  summarized  in  Table  2. 
The  order  of  the  rows  corresponds  to  the  order  of  the  calculation  steps. 
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Table  2  Solutions  to  the  normal  contact  problem  for  simple  profiles 


[0,  r<a 

3D  profile /(r)  \  r2  /  2R  rtanfl  cr" 

[oo5  r  >a 

fO,  JC<<3  K.  , 

ID  profile  may  a;  IB.  —  \x tan#  kjf,  jc 

contact  radius    2    ^  ^  J  *" 

a,  according  to  a 
(10) 


tt  tanf? 


normal  force  FN,  2aE'd  -E'RmdlrL  2  E"  -  2n    E'd  " 

according  to  (12)  '  3  7t      tailf?         „  + 1  (cKK$!s 

E'd  1 


stress  distribution 


(14) 


p<r),accordingto      ko.         (ry  \R\    f    I      I        Eq.  (16) 


1- 


displacement                     /■  \ 
u,(r)  according      — arcsin  -              Eq.  (18)              Eq.  (17) 
to(15)forr>Q  n_   


4.  Adhesive  Normal  Contact 

The  MDR  rule  for  the  mapping  of  adhesive  contacts  will  be  formulated  in  the  following. 
As  with  the  non-adhesive  contact,  the  MDR  transformed  one-dimensional  profile  is  brought 
into  contact  with  the  linearly  elastic  foundation  defined  in  Section  2.  Now,  it  will  be  assumed 
that  all  of  the  springs  in  the  contact  adhere  to  the  indenter,  then  the  contact  radius  remains 
the  same  after  a  successive  decrease  in  the  normal  force.  From  the  edge  of  the  contact 
towards  the  middle,  however,  more  and  more  springs  will  be  loaded  in  tension.  As  soon  as 
the  change  in  length  of  the  outer  springs  reaches  the  maximum  allowable  value: 


AZ(+fl)  =  AZma>):=  (21) 


there  will  be  a  state  of  indifference  between  adhesion  and  separation  (Fig.  5).  Here,  Ay  is 
the  separation  energy  of  the  contacting  bodies  per  unit  area.  This  state  corresponds 
exactly  to  the  equilibrium  state  of  the  three-dimensional  adhesive  contact  [15]. 

In  contrast  to  the  algorithm  for  non-adhesive  contact,  Equation  (10)  must  only  be 
replaced  by: 

uz  (±b)  =  -A/max(a)  =>  d  :=  g(a)  -  A/max(a)  (22) 
in  order  to  calculate  the  indentation  depth.  The  normal  force  is  given  as  before  by  Eq.  (12): 


FN  =  2E*  j(d  -  g(x))dx  =  2E*  ad      g(x)dx  =  2E*  ag(a)-jg(x)dx-aAlnmx(a) 


(23) 
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)  /)  />  /)  j)  i>  jf  i>  j)  t)  j  )  j>  / 

Fig.  5  Illustration  for  the  MDR  rule  for  an  adhesive  contact 


A'  A") 


a 


or 


FN  =2E*  $xg'(x)dx-aAlmJa) 


(24) 


0 


If  the  force  is  controlled  during  the  separation,  then  critical  value  ac  of  the  contact  radius 
at  the  moment  of  the  loss  of  stability  is  determined  by  condition  dFN/da  =  0: 


Inserting  the  critical  radius  obtained  from  this  equation  into  Eq.  (23)  results  in  the 
maximum  negative  force.  We  will  call  its  magnitude  the  adhesion  force  FA: 


The  simplest  example  is  calculating  the  adhesion  force  between  a  cylinder  with  the 
radius  a  and  an  elastic  half-space.  In  this  case,  the  integral  in  Eq.  (26)  is  equal  to  zero  and 
the  adhesion  force  is  only  given  by  the  first  term:  FA  =  2E*aAlimil  (a)  =  ^ira/lpAy  ,  which 
corresponds  to  the  result  of  Kendall  [20].  Calculations  for  other  profiles  are  just  as  simple 
and  are  summarized  in  Table  3. 

In  order  to  stress  the  generality  and  simplicity  of  the  calculation  method,  we  will  now 
conduct  the  calculation  for  a  parabolic  profile  with  a  worn  tip  (Fig.  6).  If  we  take  the 
equivalent  profile  from  Equation  (8)  into  account,  the  resulting  indentation  depth  is: 


(25) 


(26) 


0 


d  (a)  :=  g(a)  -  A/max  (a)  =  ^  si  a     h'  -  l-^ 


(27) 


Fig.  6  Qualitative  presentation  of  the  adhesive  contact  of  a  parabolic  profile  with  a  flattened  tip 
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and  from  Eq.  (12),  the  normal  force  is: 

FN{a)  =  E'  f  uz(x)dx  =  (2a2  +b2)4a  -b2  -^xaE'Ay  .  (28) 

J  3R 

-a 

This  equation  corresponds  to  the  results  from  the  three-dimensional  theory  [22]. 


Examples  for  adhesive  normal  contacts 

Inserting  the  MDR  transformed  profiles  into  Eqs.  (25)  and  (26)  and  conducting  an 
elementary  integration  provides  the  results  summarized  in  Table  3  for  the  "classical  profiles" 
of  the  cylinder  [20],  sphere  [19],  and  cone  [21],  as  well  as  the  general  power  function  profile 
[10,  11].  The  order  of  the  rows  corresponds  to  the  order  of  the  calculation  steps. 


Table  3  Solutions  for  adhesive  contacts  for  simple  profiles 
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5.  Tangential  Contact 

We  now  consider  an  axially-symmetric  indenter  that  is  initially  pressed  into  an  elastic 
half-space  with  normal  force  FN  and  subsequently  loaded  with  tangential  force  Fx  in  the  x- 
direction  (Fig.  7).  We  assume  that  Coulomb's  law  of  friction  is  valid  in  the  simplest 


Fig.  7  Qualitative  presentation  of  the  tangential  contact 
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form  in  the  contact:  as  long  as  tangential  stress  z  is  smaller  than  the  coefficient  of  friction 
jj,  times  normal  stress  p,  the  surface  is  in  a  state  of  sticking.  After  slip  sets  in,  the 
tangential  stress  remains  constant  and  equal  to  jjp: 

r(r)<jup(r)  for  stick,  (29) 

r(r)  =  /Up(r)    for  slip.  (30) 

It  is  known  that  for  a  small  tangential  force  at  the  edge  of  the  contact,  a  ring-shaped 
slip  domain  develops,  which  expands  inwards  for  increasing  force  until  the  complete  slip 
is  exhibited.  We  denote  the  inner  radius  of  the  slip  domain  (or  the  radius  of  the  stick 
domain)  with  c. 

The  MDR  is  applied  to  the  tangential  contact  as  follows.  The  modified  profile  g(x)  is 
pressed  into  the  linearly  elastic  foundation  with  force  FN  and  then  tangentially  displaced 
by  ux0>.  The  linearly  elastic  foundation  is  denoted  by  the  stiffnesses  according  to  (1)  and 
(2).  Every  spring  sticks  to  the  indenter  and  is  displaced  along  with  it  as  long  as  the 
tangential  force  AFX  =  kx  uj0>  is  smaller  than  fjAFz.  After  the  adhesion  force  is  reached, 
the  spring  begins  to  slip  and  the  force  remains  constant  and  equal  to  f/AFz.  This  rule  can 
also  be  incrementally  formulated  so  that  it  can  be  applied  for  arbitrary  loading  histories: 
for  a  small  displacement  of  indenter  of  Am./0',  we  obtain: 

Aux  (x)=Aux0),        if  | kxux  (x)|  <  jufz 

,  ,     ,  fAFAx)    .  -  •  <31) 

ux  (x)  =  +  ,  in  a  state  of  slip 

*, 

The  sign  of  the  last  equation  is  dependent  on  the  direction  of  motion  of  the  indenter.  By 
following  the  incremental  changes  in  the  position  of  the  indenter,  we  can  explicitly  determine 
the  displacement  of  all  the  springs  in  the  contact  area;  with  this,  all  tangential  forces: 

AFX  =  kxux (x)  =  G*Axux (x) .  (32) 

and  the  linear  force  density  (distributed  load): 

qx{x)  =  ^  =  G*ux{x)  (33) 
Ax 

are  also  known.  The  distribution  of  tangential  stress  z(r)  as  well  as  displacements  ux(r)  in 
the  original  three-dimensional  contact  are  determined  as  follows  [1]: 

2  r    u  (x) 

ux(r)  =  -\^U=dx.  (34) 
r(r)=  I   n  Tdx  =  J   n  (35) 


2  2 

x  — r 


If  the  indenter  is  displaced  in  one  direction  from  the  equilibrium  position,  then  radius  c  of 
the  stick  domain  is  determined  from  the  condition  that  tangential  force  kxux0)  is  equal  to  // 
times  normal  force  kzuz(c)  (Fig.  8): 
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G*uf  =  jUE\d-g(c)). 
The  tangential  displacement  is  equal  to: 

for  x  <  c 


ux  (x)  ■■ 


.CO) 


(d  —  g(x)),  for  c  <  x  <  a 


(36) 


(37) 


the  distributed  load  is: 


q(x)  ■■ 


Gw™, 


for  x  <  c 


/uE  (d  -  g(x)),    for  c  <  x  <  a 


and  the  resulting  tangential  force  is2: 

Fx  =2^qx{x)dx  =  2/uE"  c(d  —  g(c))+  [ (d  —  g(x))dx 

0  L  c 

The  normal  force  is  still  given  by  Equation  (12)  and  ratio  Fxl(/uFN)  is  given  by: 

xg'(x)dx 

17  J 


(38) 


(39) 


ag(a)-jg(x)dx 


(40) 


The  relative  displacement  Aux  =  mx(0)  -  ux(x)  of  the  surfaces  in  contact  is  obtained  by 
subtracting  ux    from  Eq.  (37): 


//////////////////////// 
Fig.  8  Equivalent  model  for  the  classical  tangential  contact  according  to  Cattaneo  and  Mindlin 


Au .  = 


0. 


for  x  <  c 


{g(c)-g(x)),  for  c<x<a 


(41) 


2  We  stress  once  more  that  all  macroscopic  values  obtained  using  the  procedure  described  above  correspond 
exactly  to  the  three-dimensional  solutions  of  Cattaneo  [4],  Mindlin  [5],  Jaeger  [6]  and  Ciavarella  [7]. 
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The  relative  displacement  in  the  original  three-dimensional  system  is  calculated  using  Eq. 
(34)  as: 


ux(r)  =  -ju 
n 


E\tg(c)-g(x) 


4 


G'){  ^r2-*2 
For  example,  for  a  conical,  we  obtain: 


dx,    for  c  <  r  <  a  . 


(42) 


ux(r)  =  ju 


tan6> 


1         •  I  c  )  n 
-1  -arcsin  -  +  — 

cj  UJ  2 


for  c<r<a.  (43) 


Examples  for  tangential  contacts 

Inserting  the  MDR  transformed  profiles  into  Eqs.  (36)  and  (40),  results  in  the 
relationships  between  radius  c  of  the  stick  domain,  ratio  Fx  /(juFN),  and  tangential 
displacement  ux0).  The  results  for  the  "classical  profiles"  of  the  sphere  [4,  5],  cone  [23], 
as  well  as  the  power  function  profile  are  summarized  in  Table  4.  The  order  of  the  rows 
corresponds  to  the  order  of  the  calculation  steps. 

Table  4  Solutions  for  the  tangential  contacts  of  simple  profiles 


3D  profile  f(r) 

ID  profile  g(x) 

«(c)  according  to 
(36) 

according  to  (40) 

Maximum  static 

displacement 
according  to  (36) 


r1  l?R 
x  j  R 
E' 


R 


E 


1- 


rtan0 


\x  tan# 


d  ctand 

2 


1- 


1- 


»d(E'/G')  fid(E'/G*} 


Hd{E'/G') 


6.  Conclusions 

In  the  present  paper,  we  have  limited  ourselves  to  the  essential  rules  and  procedures  of 
the  method  of  dimensionality  reduction.  Evidence  for  the  statements  herein  can  be  found 
in  the  works  [1]  and  [15]. 

The  possibilities  of  the  MDR  are  much  more  expansive  as  presented  in  this 
composition.  Further  successful  applications  have  been  found  for  the  rolling  contact  [24, 
25],  contacts  with  elastomers  [26],  contacts  of  rough  surfaces  [27,  28],  elastomer  friction 
[29],  thermal  effects  in  contacts  [1],  acoustic  emission  in  rough  contacts  [31],  and  wear 
[32].  Interested  readers  are  referred  to  the  cited  literature  as  well  as  the  book  [1].  In  many 
cases,  the  MDR  allows  an  analytical  solution  to  the  problem,  as  shown  in  this  work. 
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However,  it  can  also  be  easily  implemented  numerically  and  used  for  the  investigation  of 
systems  with  complex  dynamic  loadings  [33]. 

Let  us  stress  that  the  presented  form  of  the  MDR  is  only  applicable  to  contacts  with 
homogeneous  elastic  or  viscoelastic  half-spaces  and  it  does  not  take  into  account  size  effects 
[34].  However,  extensions  to  contacts  with  final  bodies  or  heterogeneous  media  are  also 
possible  [35].  The  application  of  the  MDR  to  rough  contacts  requires  a  separate  paper. 
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METODA  DIMENZIONALNE  REDUKCIJE  U  KONTAKTNOJ 
MEHANICI I  TRIBOLOGIJI:  UPUTSTVO  ZA  KORISNIKE 
I.  AKSIJALNO-SIMETRICNI  KONTAKTI 

Metoda  dimenzionalne  redukcije  (MDR)  predstavlja  metodu  za  proracun  i  simulaciju  kontakta 
elasticnih  i  viskoelasticnih  tela.  Sustinski  se  sastoji  od  dva  jednostavna  koraka:  (a)  zamene 
trodimenzionalnog  kontinuuma  jedinstveno  definisanom  jednodimenzionalnom  linearno  elasticnom  Hi 
viskoelasticom  podlogom  (Winklerova  podloga)  i  (b)  transformacije  trodimenzionalnog  profila 
kontaktnih  tela  pomocu  MDR-transformacije.  Nakon  izvrsenja  ova  dva  koraka,  kontaktni  problem 
se  moze  smatrati  resenim.  Za  aksijalno-simetricne  kontakte  samo  mali  deo  proracuna  prevazilazi 
granice  elementarnog  racuna  i  taj  deo  nece  predstavljati  prepreku  za  prakticno  orijentisanog 
inzenjera.  Kao  alternativno  resenje,  MDR  se  moze  primeniti  i  numericki,  sto  je  gotovo  trivijalno 
zbog  nezavisnosti  elemenata  podloge.  I  pored  jednostavnosti,  svi  rezultati  su  egzaktni.  Rad 
predstavlja  kratak prakticni  vodic  kroz  MDR. 

Kljucnereci:  normalni  kontakt,  tangencijalni  kontakt,  adhezija,  trenje,  delimicno  klizanje,  napon 


